Modeling human dynamics of face-to-face interaction networks 
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Face-to-face interaction networks describe social interactions in human gatherings, and are the 
substrate for processes such as epidemic spreading and gossip propagation. The bursty nature of 
human behavior characterizes many aspects of empirical data, such as the distribution of conver- 
sation lengths, of conversations per person, or of inter-conversation times. Despite several recent 
attempts, a general theoretical understanding of the global picture emerging from data is still lack- 
ing. Here we present a simple model that reproduces quantitatively most of the relevant features 
of empirical face-to- face interaction networks. The model describes agents which perform a random 
walk in a two dimensional space and are characterized by an attractiveness whose effect is to slow 
down the motion of people around them. The proposed framework sheds light on the dynamics 
of human interactions and can improve the modeling of dynamical processes taking place on the 
ensuing dynamical social networks. 
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Uncovering the patterns of human mobility [1] and so- 
cial interactions [2] is pivotal to decipher the dynamics 
and evolution of social networks [3], with wide practi- 
cal applications ranging from traffic forecasting to epi- 
demic containment. Recent technological advances have 
made possible the real-time tracking of social interac- 
tions in groups of individuals, at several temporal and 
spatial scales. This effort has produced large amounts 
of empirical data on human dynamics, concerning letter 
exchanges 1 4] , email exchanges mobile phone commu- 
nications [1], or spatial mobility [6|, among others. 

Specially noteworthy is the data on face-to-face hu- 
man interactions recorded by the SocioPatterns collab- 
oration [7] in closed gatherings of individuals such as 
schools, museums or conferences. SocioPatterns deploy- 
ments measure the proximity patterns of individuals with 
a space-time resolution of ~ 1 meter and ~ 20 seconds 
by using wearable active radio-frequency identification 
(RFID) devices. The data generated by the SocioPattern 
infrastructure show that human activity follows a bursty 
dynamics, characterized by heavy-tailed distributions for 
the duration of contacts between individuals or groups of 
individuals and for the time intervals between successive 
contacts [1, IH . These observations contradict traditional 
frameworks positing Poisson distributed processes, and 
corroborate previous evidence of bursty human dynam- 
ics based on different datasets 0, [HI . 

The bursty dynamics of human interactions has a deep 
impact on the properties of the temporally evolving net- 
works defined by the patterns of pair-wise interactions 
[iq| , as well as on the behavior of dynamical processes 
developing on top of those dynamical networks 0, [lll - 
[l6| . A better understanding of these relevant issues calls 
for new models, capable to reproduce the bursty char- 
acter of social interactions and trace back their ultimate 



origin, beyond considering their temporal evolution jlTj . 
Previous modeling efforts mostly tried to connect the ob- 
served burstiness to some kind of cognitive mechanisms 
ruling human mobility patterns, such as a reinforcement 
dynamics [l8|, cyclic closure [l9| or preferential return 
rules [ioj , or by focusing on the relation between activity 
propensity and actual interactions [l7| . 

In this Letter, we present a simple model of mobile 
agents that captures most distinctive features of the em- 
pirical data on face-to-face interactions recorded by the 
SocioPatterns collaboration. Avoiding any a priory hy- 
pothesis on human mobility and dynamics, we assume 
that agents perform a random walk in space [2l| and 
that interactions among agents are determined by spa- 
tial proximity [22| . The key ingredients of the model are 
the following: We consider that individuals have differ- 
ent degrees of social appeal or attractiveness^ due to their 
social status or the role they play in social gatherings, as 
observed in many social [23|, economic [24] and natural 
(iBj communities. The effect of this social heterogeneity 
is that interactions, as well as the random walk motion of 
the agents, are biased by the attractiveness of the peers 
they met over time. Additionally, we assume, according 
to experimental data, that not all the agents are simulta- 
neously present in system, but can jump in and out of an 
active state in which they can move and establish inter- 
actions. We will see that these simple assumptions allow 
the model to reproduce many of the distinctive properties 
of face-to-face human interaction networks. 

The model is defined as follows (see Fig. [T]): N agents 
are placed in a square box of linear size L with peri- 
odic boundary conditions, corresponding to a density 
p = N / L'^ . Each individual i is characterized by her 
attractiveness or social appeal, which represents her 
power to raise interest in the others. The attractiveness 



ai of the agents is a (quenched) variable randomly chosen 
from a prefixed distribution 77(a), and bounded in the in- 
terval ai G [0, 1). Agents perform a random walk biased 
by the attractiveness of neighboring individuals. When- 
ever an agent intercepts, within a distance smaller than 
or equal to another individual, they start to interact. 
The interaction lasts as far as distance between them is 
smaller than d. Crucially, the more attractive an agent 
j is (the largest her attractiveness a^), the more interest 
she will raise in the other agent who will slow down her 
random walk exploration accordingly. This fact is taken 
into account by a walking probability Pi{t) which takes 
the form: 



Pi{t) = 1 - niax {a^}. 



(1) 



where JVi{t) is the set of neighbors of agent i at time i.e. 
the set of agents that, at time t, are at a distance smaller 
than or equal to d from agent i. Hence, the biased ran- 
dom walk performed by the agents is defined as follows: 
At each time step each agent i performs, with proba- 
bility Pi{t)^ a step of length v along a direction given by 
a randomly chosen angle ^ G [0,27r). With the comple- 
mentary probability 1 — Pi{t)^ the agent does not move. 
Thus, according to Eq. ([T]), if an agent i is interacting 
with other agents, she will keep her position in the fol- 
lowing time step with a probability proportional to the 
appeal of his most interesting neighbor. 

Furthermore, the empirical observations of SocioPat- 
terns data show that not all the agents involved in a social 
event are actually present for its entire duration: Some 
agents leave the event before the end, some join it later 
after the beginning, and some others leave and come back 
several times. Therefore we assume that agents can be in 
an active or an inactive state. If an individual is active, 
she moves in space and interacts with the other agents; 
otherwise she simply rests without neither moving nor 
interacting. At each time step, one inactive agent i can 
become active with a probability r^, while one active and 
isolated agent j (not interacting with other agents) can 
become inactive with probability 1 — rj. The activation 
probability of the individual i thus represents her ac- 
tiveness in the social event, the largest the activity r^, 
the more likely agent i will be involved in the event. We 
choose the activation probability of the agents ran- 
domly from an uniform distribution C{r)^ bounded in 
ri G [0,1], but we have verified that the model behav- 
ior is independent of the activity distribution functional 
form (even if we consider a constant activity rate, = r 
for all agents, we obtain very similar results). 

Within this framework, each individual performs a dis- 
crete random walk in a 2D space, interrupted by interac- 
tions of various duration with peers. The model is Marko- 
vian, since agents do not have memory of the previous 
time steps. The full dynamics of the system is encoded in 
the collision probability Pc = pird'^ , the activation prob- 
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FIG. 1. (color online) Left: Blue (dark) colored agents are ac- 
tive, grey (light) agents do not move nor interact. Interacting 
agents, within a distance are connected by a link. Right: 
The interactions of the agent represented at the center of a 
disk of radius d are presented. Each individual is character- 
ized by a number representing her attractiveness. The proba- 
bility for the central individual to move is p — 1.0 — 0.6 = 0.4, 
since the attractiveness of the inactive agent is not taken into 
account. 



ability distribution C{^)^ the attractiveness distribu- 
tion r]{a). The latter can hardly be accessed empirically, 
and is likely to be in its turn the combination of different 
elements, such as prestige, status, role, etc. Moreover, in 
general attractiveness is a relational variable, the same 
individual exerting different interest on different agents. 
Avoiding any speculations on this point, we assume the 
simplest case of a uniform distribution for the attractive- 
ness [Isl . Remarkably, this simple assumption leads to a 
rich phenomenology, in agreement with empirical data. 

In the following we will contrast results obtained by 
the numerical simulation of the model against empir- 
ical results from SocioPatterns deployments in several 
different social contexts: a Lyon hospital ("hosp"), the 
Hypertext 2009 conference ("ht"), the Societe Fran- 
caise d'Hygiene Hospitaliere congress ( "sfhh" ) and a high 
school ("school"). A summary of the basic properties of 
the datasets is provided in Table |T] (see Refs. ]i, §, ^ 
for further description and details). The model has been 
simulated adopting the parameters v = d = 1^ L = 100 
and N = 200. Different values of the agent density p 



Dataset 


N 


T 


P 


(At) 


(k) 


(s) 


hosp 


84 


20338 


0.049 


2.67 


30 


1145 


ht 


113 


5093 


0.060 


2.13 


39 


366 


school 


126 


5609 


0.069 


2.61 


27 


453 


sfhh 


416 


3834 


0.075 


2.96 


54 


502 



TABLE I. Some properties of the SocioPatterns datasets un- 
der consideration: N, number of different individuals engaged 
in interactions; T, total duration of the contact sequence, in 
units of the elementary time interval to = 20 seconds; p, aver- 
age number of individuals interacting at each time step; (At), 
average duration of a contact; (k) and (s): average degree 
and average strength of the projected network, aggregated 
over the whole sequence (see main text). 
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are obtained by changing the box size L. In the initial 
conditions, agents are placed at randomly chosen posi- 
tions, and are active with probability 1/2. Numerical 
results are averaged over 10^ independent runs, each one 
of duration T up to T^ax = 2 x 10^ time steps. 

The temporal pattern of the agents' contacts is prob- 
ably the most distinctive feature of face-to-face interac- 
tion networks [sl, [9[. We therefore start by considering 
the distribution of the duration At of the contacts be- 
tween pairs of agents, P(At), and the distribution of 
gap times r between two consecutive conversations in- 
volving a common individual, P(r). The bursty dynam- 
ics of human interactions is revealed by the long-tailed 
form of these two distributions, which can be described 
in terms of a power-law function [H . Figures [2] and [3] 
show the distribution of the contacts duration P{At) and 
gap times P(r), respectively, for the various sets of em- 
pirical data along with the same distributions obtained 
by simulating the model described above with density 
p = 0.02. In the case of the contact duration distri- 
bution, numerical and experimental data match almost 
perfectly, see Fig. [2jmain). Moreover, numerical results 
are robust with respect to variations of the collision prob- 
ability Pc = 7T(Pp^ see Fig. [2jinset). It also worth high- 
lighting the crucial role played by the heterogeneity of 
attractiveness a^. In fact, assuming it constant, = a 
(and neglecting excluded volume effects between agents) 
our model can be mapped into a simple first passage time 
problem [Hj, leading to a distribution P{At) ^ (At)"^/^ 
with an exponential cut-off proportional to — a). 

The (non-local) convolution of the exponential tails in- 
duced by the heterogeneous distribution of attractiveness 
leads in our model to a power law form, with no apparent 
cut-off, and with an exponent larger than 3/2, in agree- 
ment with the result observed in the SocioPatterns data. 
Regarding the distribution of gap times, P(r), the model 
also generates a long-tailed form, which is compatible, 
although in this case not exactly equal, to the empirical 
data, see Fig. [3fmain). The behavior of the distribution 
P(r) yielded by the model is substantially independent 
of the agent density p also in this case, see Fig. [3finset). 

Sociopatterns data can be naturally analyzed also in 
terms of temporally evolving graphs [lo|, whose nodes 
are defined by the agents, and whose links represent in- 
teractions between pairs of agents. Instantaneous net- 
works are thus formed by isolated nodes and small 
groups of interacting individuals, not necessarily form- 
ing a clique. Integrating the information of these instan- 
taneous graphs over a given time window produces an 
aggregated weighted network where the weight Wij 
between nodes i and j represents the total temporal du- 
ration of the contacts between agents i and j. The weight 
distribution P{w) of the various datasets are broad [l,!!!, 
see Fig. SJ^main), showing that the heterogeneity in the 
duration of individual contacts persists even when con- 
tact durations are accumulated over longer time inter- 
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FIG. 2. (color online) Main: Distribution of the contact du- 
ration, P(At), for various datasets and for the attractiveness 
model. Inset: Same distribution for the attractiveness model 
with different density p. 




FIG. 3. (color online) Main: Distribution of the time interval 
between consecutive contacts, P(t), for various datasets and 
for the attractiveness model. Inset: Same distribution for the 
attractiveness model with different density p. 

vals. Fig. S] shows that the outcome of the model is again 
in excellent agreement with all empirical data, with the 
exception of the "hosp" database. The reason of the de- 
parture of this dataset with respect of both other dataset 
and the model can be found in its relative small size (see 
Table H]), which probably leads to inaccurate statistics for 
the topological properties. 

Face-to-face networks can be further characterized by 
looking at the correlation between the number of dif- 
ferent contacts and the temporal duration of those con- 
tacts. These correlations be estimated by measuring the 
strength Si of a node defined as Si = Wij and rep- 
resenting the cumulated time spent in interactions by 
individual i, as a function of its degree /c^, defined as the 
total different agents with which agent i has interacted. 
Fig. m^inset) shows the growth of the average strength 
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FIG. 4. (color online) Main: Weight distribution P{w) for 
various empirical datasets and for the aggregate network ob- 
tained by simulating the attractiveness model. Inset: Average 
strength of nodes of degree /c, as a function of /c, for vari- 
ous empirical datasets and for the aggregate network obtained 
by simulating the attractiveness model. 



of nodes of degree /c, as a function of k in the em- 

pirical datasets and in the aggregated network obtained 
with the attractiveness model. As one can clearly see, all 
distributions (again with the exception of the smallest 
"hosp" dataset) are well fitted by a power law function 
s{k) ~ /c" with a > 1, with good agreement between real 
data and the model results. The observed super-linear 
behavior implies that on average the nodes with high de- 
gree are likely to spend more time in each interaction 
with respect to the low-connected individuals Q. 

A final important feature of face-to-face interactions, 
also revealed in different context involving human mo- 
bility [20|, is that the tendency of an agent to interact 
with new peers decreases in time. This fact translates 
into a sub-linear temporal growth of the number of dif- 
ferent contacts of a single individuals (i.e. the aggregated 
degree ki{t))^ k{t) ~ t^, with ji < 1. Fig. [5] shows the 
evolution of k{t) versus time for several agent with fi- 
nal aggregated degree /c(T), both belonging to a single 
dataset (main) and for the different datasets (inset). The 
sub- linear behavior of k{t) is clear, with ji = 0.6 ± 0.15 
depending on the dataset. Moreover, the shapes of the 
k(t) functions can be collapsed in a single curve by ap- 
propriately rescaling the data as k{t)/k{T) as a function 
of t/T, Fig. [5f inset). Fig. [5] shows that, remarkably, the 
attractiveness model is also capable to reproduce the be- 
havior of /c(t), up to the rescaling with total T time, again 
with the exception of the small "hosp" dataset. 

In summary, in this paper we have introduced a simple 
model of mobile agents that naturally reproduces the so- 
cial context described by the Sociopatterns deployments, 
where several individuals move in a closed environment, 
such as a conference venue, and interact between them 




FIG. 5. (color online) Main: Aggregated degree k{t) versus 
time for various individuals with different final degree /c(T), 
for the "ht" dataset (symbols) and for the aggregate network 
obtained by simulating the attractiveness model (line). Inset: 
Rescaled aggregated degree k{t)/k{T) as a function of time 
t/T for various empirical datasets and for the attractiveness 
model. 



when situated within a small distance (the exchange 
range of RFID devices). The main hypotheses intro- 
duced in the model are: (i) Agents perform a biased ran- 
dom walk in two-dimensional space; (ii) their interactions 
are ruled by an heterogeneous attractiveness parameter, 
Eq. ([1]); and (iii) not all agents are simultaneously ac- 
tive in the system. Without any data-driven mechanism, 
and with no explicit assumptions on the agents' cogni- 
tion, the model defined in terms of these three elements 
is able to quantitatively capture most of the properties of 
the pattern of interactions between agents, both at level 
of the broad distributions of contact and inter-contact 
times, and at the level of the ensuing temporal network. 
Importantly, results are robust with respect to variations 
of the model parameters, i.e., the collision probability 
Pc and the activity distribution functional form, C(^)- 
We have additionally checked that results do not depend 
qualitatively on the nature of the motion rule, given by 
Eq. ([1]). Indeed, other rules for the walking probabil- 
ity, such as considering the average of the attractives of 
the neighbors, i.e. pi{t) = 1 - Y.jeJVi{t) /^^i'^)^ ^^^^ 
substantially to the same behavior produced by Eq. ([T]), 
Overall, the proposed framework represents an impor- 
tant step forward in the understanding of face-to-face 
dynamical networks. It also opens interesting directions 
for future work, including the study of dynamical pro- 
cesses taking place on face-to-face networks and possible 
extensions of the model to more general settings. 
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